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Coherent control of individual two-level sys¬ 
tems (TLSs) is at the basis of any imple¬ 
mentation of quantum information. An im¬ 
pressive level of control is now achieved us¬ 
ing nuclearEl^, vacanciesl^ and charge spinsl^. 
Manipulation of bright exciton transitions in 
semiconductor quantum dots (QDs) is less ad¬ 
vanced, principally due to the sub-nanosecond 
dephasingEl. Conversely, owing to their robust 
coupling to light, one can apply tools of nonlin¬ 
ear spectroscopy!^ to achieve all-optical com¬ 
mand. Here, we report on the coherent ma¬ 
nipulation of an exciton via multi-wave mix¬ 
ing. Specifically, we employ three resonant 
pulses driving a single InAs QD. The first 
two induce a four-wave mixing (FWM) tran¬ 
sient, which is projected onto a six-wave mix¬ 
ing (SWM) depending on the delay and area of 
the third pulse, in agreement with analytical 
predictions. Such a switch enables to demon¬ 
strate the generation of SWM on a single emit¬ 
ter and to engineer the spectro-temporal shape 
of the coherent response originating from a 
TLS. These results pave the way toward multi¬ 
pulse manipulations of solid state qubits via 
implementing the NMR-like control schemes!^ 
in the optical domain. 

An appealing strategy to realize optically controlled 
quantum networks in solid state, is to coherently 
couple distant TLSs - like bright exciton transi tions 
in QDs - via photons confined in microcavities^^^Iiil 
or propagating in waveguides^^^HUl. xhe efficient re¬ 
trieval and manipulation of coherent responses from 
single excitons is mandatory for assessing the prop¬ 
erties of such solid state qubits and demonstrating 
the coherence transfer between them. In this con¬ 
text, substantial progress has been made by the in¬ 
troduction of the heterodyne spectral interferometry 
technique!^ to measure four-wave mixing (FW M) fr om 
individual excitons in various nanostructures^^^^^^^^^i^. 
Yet, measuring on single strongly-confined excitons in 
InAs QDs remains challenging and has not been pre¬ 
viously achieved. This is due to their small dipole 
moment /r and the resulting required high resonant 
amplitudes of the three fields; £i, £ 2 , £ 3 , employed to 
drive the measured FWM polarization, being propor¬ 
tional to fj,'^£*£ 2£3 and higher order terms with the 
same phase dependence on the fields. 

Here, the retrieval and manipulation of wave-mixing 


signals from single InAs QDs is accomplished by em¬ 
bedding them in a low-Q planar semiconductor micro- 
cavity*^, as detailed in Supplementary Fig. S2. Such 
a semiconductor hetero-structure provides an intra¬ 
cavity field enhancement of ^/Q ~ 13 at the QDs posi¬ 
tion, whilst offering spectral matching of the photonic 
mode with femto-second laser pulses. This cavity- 
enhanced driving allows reducing the resonant exci¬ 
tation intensity, required to induce the FWM, by a 
factor around ~ 2200. In consequence, the re¬ 
trieval efficiency of wave-mixing signals from single 
excitons is improved by two orders of magnit ude w ith 
respect to previous experiments on bare QDs*^*^^. In 
this work, we employ such wave-mixing responses to 
realize a novel scheme for ultrafast coherent control 
of optically active TLSs. We demonstrate gating of 
their coherent emission, by converting the FWM po¬ 
larization into the SWM one. We show engineering of 
the FWM spectral response, also acting on the global 
spectral lineshape of a TLS. Details regarding the cur¬ 
rent experimental configuration are given in Supple¬ 
mentary Fig. SI. 

For the purpose of the experiment described in this 
Letter, we require an optically active TLS in a solid. 
In InAs QDs, this is the case of negatively charged 
excitons (trions)*^, where the level structure can be 
trimmed down to two levels (disregarding spin degen¬ 
eracy): the ground and excited states corresponding 
to the presence of an electron and a trion in the QD, 
respectively, as depicted in Fig.[2a. We have there¬ 
fore used a sample which is intentionally n-doped*^ 
with a Silicon (5-layer located in the cavity spacer. To 
illustrate the coherent response of trions, we present 
in Fig.[2a the spectral interferogram detected at the 
FWM frequency 2 O 2 — £li, as detailed in the Supple¬ 
mentary Material. Based on series of interferograms, 
we construct the FWM hyperspectral imaging™ as 
exemplified in Fig. [T]b. Therein, the brightest, local¬ 
ized peaks correspond primarily to the FWM gener¬ 
ated by trions. Such imaging is employed to determine 
spatial and spectral position of QDs. Also it enables a 
detailed statistical analysis on excitonic complexes, as 
exemplified for biexcitons in Supplementary Fig. S3. 

We first retrieve the required parameters to infer 
and model the optical response of single emitters in a 
solid. The three-beam FWM spectroscopy offers di¬ 
rect access to coherence and population dynamics in 
the TLS. To probe the population dynamics, we de¬ 
tect the time-integrated FWM at the UI 3 Ut 2 — 
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FIG. 1. Wave mixing spectroscopy of individual 
negative trions in InAs QDs embedded in a low- 
Q microcavity, a) Interferograms of four-wave mixing 
(FWM) of a single trion for 20,2 — tli, ti 2 = 0.5 ps 
and pulse areas ( 0 i, 62 ) = (rr/Z, tt), and six-wave mixing 
(SWM) at 2 O 3 — 2 O 2 + tti for ri 2 = 0.5 ps, r 23 = Ips, 
( 61 , 62 , Oz) = ( 7 r/ 2 , TT, tt). Integration time 20s. Het¬ 
erodyned spectral shape of £2 is shown with a purple, 
dashed line, b) Amplitude of the FWM hyperspectral 
imaging at ri 2 = 1 ps. Localized peaks correspond to 
the FWM of single QDs: the amplitude is displayed 
as brightness, the energy is indicated as the hue of the 
colour range (indicated by the bar) from -0.3 (red) to 
-1-0.3 meV (blue) with respect to the chosen center of 
hoj = 1350 meV. c) Population dynamics (brown squares) 
of the trion probed with FWM, ri 2 = 0.5 ps: exponen¬ 
tial decay (brown line) yields the population lifetime of 
Ti = (390 ± 10) ps. Coherence dynamics (black circles) of 
the same trion probed with FWM, consistent with radia- 
tively limited dephasing T 2 = 2 Ti and a residual inhomo¬ 
geneous broadening of ct = (3 ± 1) /reV. d) FWM power as 
a function of \/P 7 = |£^i| and of £1 pulse area 81 showing 
Rabi oscillations, 62 = 2 / 37 r, P 2 = |£ 2 p = 157W/cm^, 
ri 2 = 0.5ps. Pulse areas of {81,62) = ( 7 r/ 2 , tt ), corre¬ 
sponding to (Pi, P 2 ) = (72, 235) W/cm^, are determined 
by the first maximum of the FWM versus 81 {62, not 
shown). Predicted sin^(0i) exp(—^^i) is depicted by the 
orange dashed-dotted line, where ^ is the damping con¬ 
stant. 


frequency, while varying the delay T 23 between £2 and 
£3 and fixing ri 2 = 0.5 ps. The resulting evolution 
of the FWM amplitude is shown in F ig.lTI c and ht- 
ted by an exponential decay, yielding^S^ the life¬ 
time of Ti = (390 ± 10) ps. Conversely, detect¬ 
ing the FWM at the frequency 2 Q 2 — fli, reflects 
the coherence in the TLS. Its evolution is governed 
by the delay ri 2 between £i and £ 2 , as displayed in 
Fig.gc. The lack of FWM at negative delays and the 
absence of fine-structure beating^ confirm the tri- 


onic nature of the investigated transition. Our data 
are well described by the product of an exponential 
and a Gaussian decay^^. From the former, we infer 
the dephasing time T 2 and the related homogenous 
broadening 72 = 2 h/T 2 (FWHM). The latter yields 
the inhomogeneous broadening a due to the residual 
spectral wandering, occurring during the integration 
time. Note that even on the single transition level 
such wandering creates a photon echo, manif ested here 
by a Gaussian decay of the coherence^^^^^. The co¬ 
herence dynamics can be fitted using T2=2Ti=780ps 
and cr = (3 ± 1) /reV. We thus find that the coher¬ 
ence of trions in these QDs approaches the radiative 
limit. This is supported by observation of their first 
order reflectance, as shown in Supplementary Fig. S5. 
More examples of such transitions close to the radia¬ 
tive limit are given in Supplementary Fig. S4. 




FIG. 2. FWM/SWM switching, a) Pulse sequence 
used in the experiment, highlighting rationale of the 
FWM/SWM switching. The evolution of the trion’s Bloch 
vector (red) during FWM, switching and SWM is also de¬ 
picted (ground state on the top of the sphere). The FWM 
transient (blue line) is created by £1 and £2 fields arriving 
at t = —T 12 and t — 0, respectively and it freely evolves 
during T23. £3, arriving at t = r23, converts the FWM into 
the SWM (green line). The FWM/SWM conversion effi¬ 
ciency reaches unity for 63 = tt . FWM (b) and SWM (c) 
amplitudes as a function of real time t and delay r 23 for 
63 = O.Stt. FWM (SWM) signals are present above (be¬ 
low) the diagonal giving the arrival of £3. Amplitude on 
a linear colour scale, as given. Decreased amplitudes for 
long times is due to the limited resolution of the spectrom¬ 
eter. d) Demonstration of the FWM suppression (blue) 
and the SWM build up (green trace) at the arrival of £3 for 
t = T 23 = 22 ps. The noise floor is given by a gray dashed 
line. Note the coexistence of both signals for t > 22 ps and 
corresponding residual signal in Fig.j^c, owing to 83 7 ^ tt . 


The coherent control experiment described be¬ 
low requires a calibration of the pulse areas 6 i = 
f dt fj,l£i(t)l/h, which are proportional to the square 
root of the pulse intensities P^. To illustrate this 
calibration, we present in Fig. [T]d the FW M power 
as a function of z/Pi and 9i. As expected F° l ^^ l ^^ l, 
the FWM undergoes Rabi oscillations with increasing 
9i. In addition, we observe a 0i-dependent damp- 
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ing, which is attributed to dissipative coupling with 
acoustic phonons^. In the three-beam heterodyne 
spectral interferometry technique, the additional de¬ 
gree of freedom provided by the time delay T23 al¬ 
lows to transform the FWM into SWM at the defined 
time after creating the FWM. Such SWM signals and 
their interferences with c oexist ing FWM have been 
studied in atomic physics^^SlSai Higher order wave¬ 
mixing has been also investigated in the condensed 
matter physics to explore many particle correlations 
in quantum wells^SSHSU qq recently 

inferring their non-Markovian dynamics!^. In the 
present study, we demonstrate for the first time gen¬ 
eration of SWM on single emitters. Here, SWM de¬ 
tection serves as a tool to implement coherent control 
upon the FWM transient of a TLS. We use £3 pulse 
to project the FWM into the SWM signal oscillat¬ 
ing at the 2173 — 2112 + frequency. The spectral 
interferogram measured at this frequency is shown 
in Fig.[^a, for the applied pulse sequence given in 
Fig.j^a. The FWM/SWM swapping is controlled by 
adjusting the pulse area 63 and r23, as derived in the 
analytical model presented in the Supplementary Ma¬ 
terial: FWM field £f and SWM field Eg read: 


\£f{t)\ oc sin01 sin^ ^ ( 0(t) -sin^ - T23) ) x 


2 "3 , 


a-t'12 


( 1 ) 


|£s(i)l «: sin0i sin^ y sin^ Y0(t - T 23 )e (2) 


where 0 is the Heaviside function, and t = 0 is 
defined by the arrival of £2- 

Conversion of FWM into SWM and their co¬ 
existence is experimentally demonstrated in Fig.[^ 
Therein, we present the FWM and SWM field am¬ 
plitude transients for increasing r23 at 63 = O.Stt. In 
such maps, the arrival of £3 at f = r23, generating the 
conversion, defines the diagonal. From Eqs. Q and 
(i). we expect that the FWM is present in the upper- 
side of the diagonal (t < T23), and has been converted 
to SWM in the lower-side {t > T23), as indeed mea¬ 
sured. This temporal gating of both signals is the key 
to manipulate the spectral distribution of the FWM 
from a TLS. 

We note that the conversion efficiency reaches unity 
for the 63 = TT yielding, at the arrival of £3, a complete 
suppression of the FWM. In Fig.[^ we employ specific 
pulse areas (0i, ^2) = (f, ti") in order to drive, between 
£2 and £3, the maximum polarization to the FWM, 
consistent with Eq. 0. In Fig.|^a we present EWM 
transients measured at 03 = tt for various r23. Their 
spectral amplitudes (also see Supplementary Eig. S6) 
are given in Fig. ib. The spectra show a substan¬ 
tial broadening and ringing due to the step-like sup¬ 
pression in time, leading to a sine cardinal like shape. 
FWM at 03 = 0 (topmost) is given for comparison. 
The effect, also clearly visible in Fig.j^c, is observable 
from r23 = 110 ps to T23 = 10 ps where the FWHM 
of the main peak increases from 38 /xeV (given by the 
spectral resolution) up to 700 /xeV. This represents an 
imposed broadening by two orders of magnitude with 
respect to 72. With further decrease of r23 such spec¬ 
tral broadening can reach 10 meV range, as is only 
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FIG. 3. Manipulation of the coherent response of a 
single emitter via wave-mixing switching, a) FWM 

transients generated for ri 2 = 0.2 ps, 03 = tt and r 23 as 
labeled. A complete suppression of the FWM is observed 
after the arrival of £3 (t > r23); the noise level is shown by 
the dashed yellow line. The free evolution of the FWM for 
03 = 0 is shown in the topmost panel, b) FWM spectral 
amplitudes of the data in a), showing a substantial spec¬ 
tral broadening of the FWM due to the step-like suppres¬ 
sion of the signal in time domain. Analytical predictions 
using Eq. ( |14[ ) are shown as red lines, c) FWM spectral 
amplitudes of the transients given in Fig.|^b, showing the 
evolution of the FWM lineshape with increasing T23 for 
03 = O.Stt. Logarithmic colour scale, as shown by the ver¬ 
tical bar. 


bounded by the duration of £3. Due to the limited 
amount of the emitted FWM for T23 — 0 , its obser¬ 
vation requires a large signal-to-noise ratio. Fourier 
transforming Eq. Q with respect to the real time t, 
yields the spectrally resolved FWM amplitude: 


Q 

|£/(w)| ocsin0i sin^X 

/ I -h Sin^ ^g- 2 T 2372 _ 2 gin^ ^g-r 2372 gQg AuJT 23 

V 27r(7| -h Aw2) 


1/2 


( 3 ) 


where Aw = w — ujeg ■ With the particular pulse ar¬ 
eas (01, 02 , O3) = (|, TT, tt), the spectrum takes the 

simple form g- 72 (Ti 2 +T 23 ) 

which reproduces the observed features quantitatively, 
as shown by red traces in Fig. ib. The spectral broad¬ 
ening of the central peak scales as characteristic 
of the sine cardinal lineshape. 

For intermediate values of 63, the FWM is only 
partially converted to SWM. Therefore both non- 
linearities coexist, as can be noted in Fig. id. This is 
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FIG. 4. Manipulation of the FWM with the area of 
the control pulse £ 3 . a) Measured evolution of the FWM 
transient for 6 ^ of O.drr (blue dashed), O.Ttt (green dotted) 
and TT (orange solid); for T 12 = 0 . 2 ps and r 23 = 20 ps. 
An increasing suppression of the FWM for t > r 23 with 
increasing 62 , is observed, b) FWM spectral amplitudes of 
the transients shown in a), with the corresponding theoret¬ 
ical prediction according to Eq. (3). c) Map of the FWM 
spectral amplitude as a function of 62 demonstrating a 
gradual suppression of the spectrally narrow component 
of the FWM. Delays as in a), d) Theoretical simulation of 
c). Logarithmic color scale, as shown with the colour bar. 


investigated in Fig.|^ where the FWM transients (a) 
and spectra (b) for increasing 6^ and hxed r23 = 20 ps 
are presented. Owing to the increasing FWM sup¬ 
pression for t > 20 ps, the amplitude of the spec¬ 
trally narrow component is reduced (see Supplemen¬ 
tary Fig. S 7 ). In parallel, the broad pedestal devel¬ 
ops. Hence, non-natural lineshapes can be designed 
by tuning 63 and r23. Eq. ( 14 ), derived in the Supple¬ 
mentary Material, reproduces the experimental data 
shown in Fig.j^b and c, without any free parameters, 
apart from the absolute common scaling. 


The ability to manipulate the spectral width of a co¬ 
herent response from the TLS provides a novel degree 
of freedom in solid state quantum optics. We empha¬ 
size that the phase shift induced by the FWM/SWM 
conversion modihes the coherent oscillation of the 


dipole: the global lineshape of the TLS polarization is 
altered, as highlighted in Supplementary Fig. S8. The 
latter is strikingly manifested when the £3 is hetero¬ 
dyned at 2fl2 — Hi, leading to a stationary phase shift 
between the two non-linearities generated by £2 and 
£3, respectively. Therefore by adjusting this phase 
shift with T23 (see Eqs. ( 11 ) and ( 12 ) in the Supple¬ 
mentary Material), we can perform controlled phase 
rotations of the TLS dipole. After arrival of £3, the 
phase of the TLS emission could be inverted or en¬ 
tirely frozen. Such on-demand blockade of emission 
represents a fundamental step in optical control of sin¬ 
gle TLS systems. As an example of application, let us 
note that tuning the spectral shape of a TLS could 
be employed to optimize injection and storing of sin¬ 
gle photons in optical resonators. Furthermore, max¬ 
imizing the spectral overlap within a pair of distant 
TLSs promotes their coherent coupling via propagat¬ 
ing photonsl^, which is a prerequisite to realize opti¬ 
cally controlled quantum networks in a solid. Multi¬ 
wave mixing could be also used to selectively address 
the coherence dynamics and transitions in higher man¬ 
ifolds of the Jaynes-Cummings ladder of nano photo nic 
devices operating in a strong-coupling regime^^^. 

By employing a low-Q semiconductor microcavity, 
we improved the retrieval sensitivity of coherent re¬ 
sponses from individual emitters, enabling wave mix¬ 
ing spectroscopy of individual excitons in strongly- 
confined InAs QDs. Using three beam configuration, 
we inferred a new scheme for controlling coherent evo¬ 
lution of a TLS via converting FWM to SWM. We 
have demonstrated that, via temporal gating of the 
FWM, we can control its lineshape by varying the 
amplitude and the arrival time of the gate pulse. The 
demonstration of SWM on individual emitters paves 
the way towards investigations of their non-Markovian 
dynamics^ and to study higher order correlations^^ 
involving more than two excitons within a QD. Ex¬ 
tending our proof-of-principle protocol toward control 
via multi-wave mixing is conceptually straightforward: 
conversion between wave mixing processes is deter¬ 
mined by areas and respective delays of driving pulses. 
It is also technically feasible via multiplexed digital 
heterodyning. An alluring perspective is to perform 
multi-wave mixing on radiatively coupled pairs of dis¬ 
tant excitons^^, achieving non-local quantum control 
in a solid. 
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SUPPLEMENTARY MATERIAL 
Multi-Wave Coherent Control of a Solid State Single Emitter 

F. Fras, Q. Mermillod, G. Nogues, C. Hoarau, C. Schneider, 

M. Kamp, S. Hofling, W. Langbein, and J. Kasprzak 

I. THEORY: DERIVATION OF THE FWM AND SWM RESPONSES 

In this part we derive the four-wave mixing (FWM) and six-wave mixing (SWM) responses, as experimentally 
studied in the main part of the manuscript. We consider a two level system (TLS) formed by an electron in 
the ground state |0) and a trion in the excited state |1). In the model developed below the coupling between 
the optical fields and the TLS is treated exactly. The optical response of the TLS is expressed at all orders of 
the exciting field amplitudes (not limited to the third and five orders, which are usually considered in the low 
field intensity regime). The relevant semiclassical Hamiltonian of the considered system in interaction with the 
three exciting fields is given hy H = Hq + Hl, 

Ho=hWeg\l){ll 

3 

Hl = e*G=*(*-.)+0,)|o)(i| + C.C., (4) 

where ujeg refers to the electron-trion optical transition frequency, fj, is the electric dipole of the QD, Si{t) 
determines the effective electric field amplitude of the pulse i at the QD position, represents the arrival time 
of the pulse i at the QD position and (f>i refers to the additional phase shift induced by optical heterodyning. 
The excitation pulses are assumed to be delta-function pulses, i.e., Si{t) = SiS{t — Ti). In addition to the 
evolution governed by the Hamiltonian H, the system undergoes dissipative dynamics due to the coupling with 
its environment. Therefore the state of the system p(t) is described in the density matrix formalism in the basis 
{| 0 ),| 1 )}. 

In order to relate the different multi-wave responses we derive the state of the system p(t). We assume 
that the pulse duration is much shorter than the relevant time scales of the system dynamics. Therefore the 
interaction between the system and its dissipative environment is neglected during the coupling with the optical 
excitation fields. Under this assumption, the modification of the state p{t) induced by the pulses is obtained 
analytically according to the unitary transformation U{9): 

p(r+) = C/(d.)p(r-)C/(0Q', (5) 

U{e,) = cos |[| 0 )( 0 | + | 1 )( 1 |] - zsin |[| 0 )(l|e^('^—-■) + c.c.], ( 6 ) 

where denotes the time instant just before (after) the pulse i, and 9i = dtp,Si{t)/h is the area 

of the pulse i. During the delay time between the pulses the open two-level system relaxes its density and 
polarization. 

In the following the dissipation is described in the Markov limit, the subsequent dynamics can be then given 
by a master equation in the Lindblad form 

P = -^[Ho,p]+C{p) (7) 

where C is the superoperator of dissipation. It takes into account the dissipation of energy (mainly induced 
by spontaneous emission in our experimental case) and pure dephasing processes, which are characterized 
respectively by the rates 71 and 7 ^. Given the minute values of fluctuations induced by the spectral wandering 
(as measured in Fig. 1. c), the inhomogeneous dephasing is disregarded here. Therefore C reads as: 

£{p) = 7 i[CpC' - ^ [dc,p}^] + \id[DpD - ^ {D\p]^] ( 8 ) 

with C = |1)(0| and D = |I)(I| — |0)(0|. 

By solving this equation we derive the complete optical polarization dynamics, which is proportional to the 
off-diagonal terms of the density matrix. We deduce therefore the field £qd radiated by the quantum dot 
structure: 

^QD^t) = A (1| p(t) |0) + c.c. with A = ipQCWegPK (9) 

where k represents a coefficient of proportionality including phase shift and absorption due to the propagation 
inside the structure. 
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The coherent optical response of the QD can be developed in a sum of contributions £m referring to particular 
nonlinear processes. The emitted fields £n are identified by their phase factors where the phase (j)n results 
from a linear combination of the exciting pulse phases : ipn = ^i4>i n = {oi, 02, 03} G 1 ?. 

The principle of the heterodyne spectral interferometry experiment is based on the capacity to select a 
specific En via corresponding phase 4 >ri- To do so, the phase relation between the exciting pulses is tailored by 
acousto-optics modulators operating at radio-frequencies fli. Therefore successive exciting pulses {Si} belonging 
to optical train i, characterized by a repetition rate are phase shifted by a factor EliTr- By repeating 
the experiment in time and combining optical heterodyning with the spectral interferometry we retrieve the 
amplitude and phase of the desired non-linear component £n out from the total electric field radiated by the 
QD dipole®. In the present study, the nonlinear signals of interest are the following: 

• the non-degenerate FWM E^/ beating at the heterodyne frequency = ps — 4>2 + 4 >i 

• the degenerate FWM Ef beating at the heterodyne frequency </>/ = 2^2 — Pi 

• the degenerate SWM Eg beating at the heterodyne frequency pg = ‘^Ps — ‘^p2 + Pi 

Considering the TLS in the ground state at the time , the dynamics of the optical fields are given by: 

£,,{t) = ^ sin(0i) sin(02) sin(d3)0(t - (10) 

’’ 2 

^ ( 11 ) 

£g[t) = -Asin(0i)sin2(y)sin2(^)0(t-r23)e-‘''^e*(“®®(‘-""^^+"i=)+^®\ (12) 

where 0 (t) is the Heaviside function, 72 = 7i/2 -I- 7^ and (ti, T2, T3) = (—7-12, 0 , r23) with = Tj — Ti and 
ti2 > 0 and r23 > 0 . Once Fourier transformed with respect to the time t, the spectrally resolved wave mixing 
field amplitudes read as : 


Ef{t) = Hsin(6»i)sin2(^) (^0(t) - sin2(^)0(t - T23) 


A g — '7’l272—■7*2371 


Ef{uj) = Asin(6»i)sin { — ) 


2 , 6 * 2 . (l-sin2(^)e^23(*^“-T'")) 


Q„ g-(Tl2-|-T23)72 




( 13 ) 

( 14 ) 

( 15 ) 


where Aw = ujgg — uj and the phases pn are omitted. In the equations ( | 10 [ ) and ( 13 ), one can note that the 
amplitude of the non degenerate FWM signal £f> is sensitive to both the relaxation of the polarization and 


density of the TLS. In particular, the expression ( 13 ) in the spectral domain highlights that probing the signal 


Ej!' as a function of the delay times ri2 and T23 allows to experimentally determine the characteristic rates 71 
and 72, respectively (see Fig. 1 in the main manuscript and Supplementary Fig. S 4 ). 

The relations 0 and ( | 12 [ ) are describing the core of the experiments presented in the main article. They 
explicitly illustrate the dynamics of the TLS coherences induced by the third pulse, where the FWM response 
is projected to the SWM response. This interplay between the two nonlinear processes is controlled by the area 
6*3 of the third pulse and the time delay T23 between the last two pulses. Conversely, such degree of freedom 
permits designing the spectro-temporal shape of the FWM component, as shown in the equations (11 14 ). 

We point out that the switching between FWM and SWM triggered by £3, physically modifies the coherent 
oscillation of the dipole and directly acts on the total lineshape of its emission. This is exemplified in Supplemen¬ 
tary Figure Fig. S8 (see page 9 ). Therein, we present simulated spectral lineshape of the complete QD optical 
polarization after the arrival of the second pulse £2 for ( 0 i, 02, Os) = (f, tt, tt). We observe that the spectral 
lineshape of the dipole is significantly altered (for example for r23 = 390 ps) compared to the unperturbed case 
{63 = 0 ). The spectral broadening is generated by the spectral interferences between the FWM (for t < T23) and 
SWM (for t > T23) components, which averaged over many repetitions of the heterodyne experiment produce 
such spectrum. 


II. EXPERIMENTAL DETAILS 

Our experiments are based on an extended version of the heterodyne spectral interferometry technique®. 
We use a pulse train spectrally centered at hujo = 1 . 35 el 4 , generated by a Ti:Saphire laser (Spectra-Physics, 
Tsunami Femto). The pulses are spectrally shaped and chirp-corrected by a passive pulse-shaper (G - grating, 
LS - lens, SM - spherical mirror). They are then spectrally up-shifted by three acousto-optic modulators 
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(AOMs) operating at ^^ 1 , 2,3 = (80, 79, 80.77) MHz. Thus created three driving fields £ 1 , 2.3 - marked as red, 
blue and green, respectively - acquire relative time delay ri 3 (positive for Si leading) and ri 2 (positive for 
Si leading) by passing through two mechanical delay lines. Their polarization can be controlled by a set of 
half- and quarter-plates, but in this work we employ co-linear, horizontal (in the lab-frame) polarization for 
all beams. The driving helds are recombined into the same spatial mode (yellow) and focussed on the sample 
surface via the microscope objective (MO - Olympus, LCPLN50XIR/0.65, BS beam sampler, FF far held) 
installed on a X-Y-Z close-loop piezo-stage, (XY-range of 250/im, Z- range of 450 ^m). The excitation pulses 
have duration of 350 ± 50 fs. The pulse chirp, mainly induced by a slab of Te02 crystal in the AOMs and in the 
MO, is compensated by measuring the non-resonant four-wave mixing produced at pulse overlap in an auxiliary 
GaAs sample and minimizing its duration with respect to ri 2 . The microcavity sample is installed inside a 
custom, continuous flow helium cryostat, containing motorized, closed-loop XY translation stages (range of 
10mm, accuracy of 50nm). Experiments reported here were performed at 5.1 ± 0.1 K. The far field is imaged 
into the AOM, where - using a Bragg condition - the reflected signal is mixed with the reference field Sr- 
The mixing AOM is driven at the frequency Od = aOi -|- 5 O 2 + CO 3 , where (a, b, c) e (—3, —2, —1, 0, 1, 2, 3) 
and a-|-b-|-c=l. These heterodyne frequencies are produced by the home-made three-channel, analogue mixer, 
with a spectral purity of 40 dB. They are subsequently regenerated by a lock-in amplifier attaining spectral 
purity of 100 dB. Mixed beams Wa and Wb are imaged into the entrance slit of the imaging spectrometer of 
750mm focal length. Spectral interferograms are retrieved by subtracting and 7r-phase flipping Wa and Wb, 
as described in Ref. [8]. We use a CCD camera from Princeton Instruments (PIXIS:400BR eXcelon), offering 
enhanced full well capacity, readout rate and quantum efficiency of 83% at 915nm. The wave mixing signal 
is retrieved in amplitude and phase by performing spectral interferometry and by adjusting Sr to arrive prior 
to the signal. We note that we have also tested optical heterodyning by performing digital mixing with an 
arbitrary wave-form generator (Tektronix, sampling rate 1.2C/s) and retrieved the FWM from single quantum 
dots with the same signal-to-noise ratio, as with the analogue mixing. This paves the way towards heterodyne 
n-wave mixing experiments, enabling quantum control of high order nonlinear processes of individual optical 
transitions in solids. 




CCD 


Supplementary Figure SI. Schema of the experimental setup. 
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III. AUXILIARY RESULTS 


The sample has been grown by the molecular beam epitaxy. A layer of annealed and capped InAs QDs 
(density 2 x 10® cm“^) is placed in a center of a GaAs spacer. A J-doping layer with Si (density 1.8 x 10^® cm“®) 
is present 10 nm below the QD layer, inducing intentional negative doping of QDs. The spacer is sandwiched 
between two Bragg mirrors. The bottom (top) mirror contains 24 (5) pairs of A/4 Alo.gGaAs/GaAs stacks, 
forming an asymmetric A-cavity with a quality factor of Q ~ 170. An example of the micro-photoluminescence 
hyperspectral imaging frame performed on the studied microcavity is shown in Fig. |S2| a. The sample is excited 
in continuous wave at 1.734 eV with 0.3/iW at the sample surface. Blue peaks correspond to the emission 
of individual InAs quantum dots. Count rates at the peak as high as 10® per seconds are routinely detected 
using a grating of 1800 grooves/mm with Low Noise and gain 2 settings of the Pixis camera. Linear color scale 
from 0 to 10® counts per seconds. In Fig. [S^b we present a typical micro-photoluminescence spectrum (black) 
and photoluminescence spectrum averaged over 400 /rm^ (gray). The reflectivity spectrum measured with the 
white light (blue) shows the cavity mode. The mode is broadened at the high energy side due to high NA of 
the microscope objective and the in-plane dispersion of the cavity. In Fig.j^c we present topography of the 
sample surface, revealed with the atomic force microscopy. Typical topography images over the sample area 
of 20 /rm x 20 /rm (left) and of 5 /im x 5 /im are shown, revealing mosaic-like morphology with characteristic 
oval photonic defects of typically 2 /im x 1 /im in-plane size and of 20 nm hight occurring with a spatial density 
0.1/im“^. Such defects offer an enhanced in- and out-coupling efficiency of around 40%, as recently shown in 
Ref. [20]. This, combined with an intra-cavity amplification, results in an enhanced excitation and detection of 
the coherent nonlinear responses of QDs situated close to the maximum of the localized optical mode induced by 
a defect. We systematically measure shorter Ti times (150-350ps), weaker driving powers and stronger FWM 
from QDs residing in the defects, with respect to the dots placed outside of defects. We speculate that these 
observations are due to the change of the photon density of states by the surface modulation defects, creating 
a varying coupling efficiency and Purcell effect. Further research is required to make a clear assignment and to 
elucidate the influence of these photonics defects on the coherent response of QDs. 

Typical frames of the FWM hyperspectral imaging at ri 2 = -flps and oj as labeled are shown in Fig.[S3|a-d. 
Linear colour scale from black to white, as shown by the vertical bar. Strongest peaks primarily correspond 
the FWM of negative trions. In order to obtain spectrally-resolved FWM image, as exemplified in Fig.j^b, 
the following procedure is applied. For each spatial position, the spectrally-resolved FWM amplitudes are 
fitted with a Lorentzian function within a chosen energy region (specifically, from 1349.7 meV to 1350.3 meV 
in Fig.[^b). The fitted center position yields the energy, while the fitted Lorentzian amplitude is proportional 
to the FWM amplitude. The retrieved parameters are rearranged into a pair of two-dimensional spatial maps; 
encoding spectral position and amplitude, respectively. They are hnally merged into a single, spatial map, 
where the amplitude is represented as brightness and the energy as a hue level. 

In Fig. S3e and f we provide the result of the FWM spectral correlation analysis AC {Sou) performed on the 


hyperspectral imaging FWM data, as defined in Ref. [38]. For ti 2 = -1-1 ps, AC{ 6 u!) reveals two peaks at Suj of 
2.2 meV and 3.8 meV attributed to neutral exciton and biexciton transitions, present in around 10% of quantum 
dots. Due to the weak strength of two particle states, specifically negatively charged biexcitons, the FWM of 
trions is vanishing at negative delays ti 2 < 0. Therefore for ti 2 = —2ps, AC{ 6 oj) shows only one peak at 
6 uj = 3.6 meV, attributed to biexciton transitions of neutral excitons. Lower signal to noise ratio is due to 
smaller statistics available. Wave mixing spectroscopy on such exciton-biexciton-trion systems will be presented 
in a forthcoming publication. 

In Fig. |S4| we provide two supplementary FWM measurements performed on two different negative trions, 
demonstrating their raditively limited dephasing (T 2 = 2 Ti). The population dynamics measured with FWM is 
shown by brown squares. Fitted mono-exponential decay (brown line) yields population lifetime of Ti=(400 ± 
7) ps (a) and Ti=(385 ± 25) ps (b). The coherence dynamics measured by FWM is depicted with black circles 
in (c) and (d), respectively. The data are fitted by the model presented in Ref. [19] (namely Eq. (2)), using 
radiatively limited dephasing and inhomogeneous broadening due to residual spectral wandering a = 7 /teV (c) 
and a = 3/xeV (d). Note the initial rise in the coherence dynamics observed in (c), due to the resolved photon 
echo formation. This figure is supplementary to Fig. 1. Pulse sequences employed to measure coherence (right) 
and population (left) dynamics are also depicted. The reference pulse £r arrives a few ps prior to t = 0. 

Being radiatively limited, such trions are expected to be observable in the first-order absorption experiment. 
In Fig. S5 we compare the FWM spectrum (blue line) with a direct reflectance spectrum of the £3 beam (black 
line) measured at the same spot. The FWM reveals dominating transition at 1367.1 meV. The reflectance 
displays a dip at the same spectral position. It corresponds to the Hrst-order absorption of the QD, with 
a contrast of a few per-cent resulting from the convolution of the homogenous linewidth with the spectral 
resolution. To further confirm that the reflectance dip indeed corresponds to the linear absorption of a trion 
in a single QD, we have performed the pump-probe experiment. The result is given in the inset, where the £2 
reflectance is shown while increasing the pulse area of the £1 beam from di = 0 (black trace) to 61 = tt (orange 
trace) for T 12 = 2ps (£1 arriving first). As expected, the £2 absorption is gradually decreased with increasing 
9i and is entirely saturated for the 9i = tt. 

In Fig.we present the amplitude of the FWM transient for T 12 = 0.2 ps, demonstrating a complete sup- 
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Supplementary Figure S2. Sample characterization: micro-photoluminescence, reflectivity and surface to¬ 
pography. 


pression of the FWM at the arrival of the £ 3 , which converts the FWM polarization into the SWM one. The 
pulse area of £3 is set to tt. Right: Measured FWM spectrum, showing the corresponding change of the FWM 
spectral lineshape depending on T 23 . The data are supplementary to Fig. 3. In Fig.[^ the FWM is generated 
only within a time window of r 23 = 6.5 ps prior to the arrival of £ 3 . As a result, with increasing 9^, the FWM 
gets virtually suppressed. The data are supplementary to Fig. 4. 
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Supplementary Figure S3. FWM hyperspectral imaging and correlation analysis. 


<^1 ^2 ^3 



0 200 400 600 0 200 400 


delay ^^ 23 (ps) 


delays12 (ps) 


Supplementary Figure S4. Radiatively limited dephasing measured on individual negative trions. Applied 
pulse sequences are depicted at the top. 
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Supplementary Figure S5. Signatures of single trions observed in the direct reflectance and the corresponding 
pump-probe spectra measured in reflectance. 
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Supplementary Figure S6. Control of the FWM transient with the delay T23, between £2 of £3. 



Supplementary Figure S7. Control of the FWM with the pulse area of £3. 
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Supplementary Figure S 8 . Modification of the global lineshape of a TLS via FWM/SWM switching. Spectrum 
of the complete dipole polarization (black solid line) after the arrival of the second pulse £2 for {Q\, 62 , ds) = (f, tt, tt), 
T 2 = 780 ps and T 23 = 390 ps. For such value of r 23 , the integrated amplitudes of the FWM and SWM components are 
almost equivalent, which maximizes the spectral interferences between them. FWM (red dashed line) and SWM (blue 
dotted line) spectra are presented as well. 







